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Abstract 

We present a description of computational program (written in Maple) for calcu- 
lation of fluxbrane polynomials corresponding to classical simple Lie algebras. These 
polynomials define certain special solutions to open Toda chain equations. 



1 Introduction 

In this paper we deal with a set of equations 

d ( z d 



h s \=p s \\h sI a -\ (i.i) 



dz \H S dz 

with the following boundary conditions imposed: 

H s (+0) = 1, (1.2) 

8 = 1, ...,r . Here the functions H s (z) > are defined on the interval (0, +oo) , P s > for all 
s and (A ss i) is the Cartan matrix for some finite dimensional simple Lie algebra Q of rank r 
( A ss = 2 for all s ) . 

The functions H s {z) > appear as moduli functions of generalized fluxbrane solutions ob- 
tained in [1] . Parameters P s are proportional to brane charge density squared Q 2 S and z = p 2 , 
where p is a radial parameter. The boundary condition (II. 2D guarantees the absence of singularity 
(in the metric) for p = +0. For fluxbrane solutions see p], [3]-[TT] and references therein. (The 
more general classes of solutions were described in |12} I13|). The simplest "fluxbrane" solution is a 
well-known Melvin solution [2] describing the gravitational field of a flux tube. The Melvin solution 
corresponds to the Lie algebra A± = sl(2) of rank 1 . 

It was conjectured in [1] that eqs. (jl.ip . (jl.2p have polynomial solutions 

n s 

H s (z) = l + Y, P s k)zk > (1-3) 

k=l 

where are constants, k = 1,. . . ,n s . Here pj" 8 -* ^ and 



2^2 A ss '. (1.4) 
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s = 1, r , where (A ss ') = (A, s /) _1 . Integers n s are components of the so-called twice dual Weyl 
vector in the basis of simple roots [J3]. It was pointed in pQ that the conjecture on polynomial 
structure of H s (suggested originally for semisimple Lie algebras) may be verified for A n and C n 
Lie algebras along a line as it was done for black-brane polynomials from [15] (see also p3]). In pQ 
certain examples of fluxbrane solutions corresponding to Lie algebras Ax © . . . © A\ and A^ were 
presented. 

The substitution of (jl.3p into (jl.ip gives an infinite chain of relations on parameters P^ and 
P s . The first relation in this chain 

P S = PW, (1.5) 

s = 1, ...,r , corresponds to z°-term in the decomposition of (jl.ip . 

Special solutions. We note that for a special choice of P s parameters: P s = n s P , P > , 
the polynomials have the following simple form [1] 

H s (z) = (l + Pz) n °, (1.6) 

s = l,...,r. This relation may be considered as nice tool for verification of general solutions 
obtained by either analytical or computer calculations. 

Remark: open Toda chains. It should be also noted that a set of polynomials H s define a 
special solution to the open Toda chain equations [161 [T71 [TBI [T9] corresponding to the Lie algebra 
G 

f^ = -B s exp(j2A ss ><f'), (1.7) 

s'=l 

where B s = 4P S , 

H s = exp(-q s (u) - n s u), (1.8) 

s = 1, r and z = e~ 2u . 

In this paper we suggest a computational program for calculations of polynomials corresponding 
to classical series of simple Lie algebras. 

2 Cartan matrices for classical simple Lie algebras 

Here we list, for convenience, the Cartan matrices for all classical simple Lie algebras and inverse 
Cartan matrices as well. 

In summary [14], there are four classical infinite series of simple Lie algebras, which are denoted 

by 



A r (r > 1), B r (r > 3), C r (r > 2), D r (r > 4). (2.1) 

In all cases the subscript denotes the rank of the algebra. The algebras in the infinite series 
of simple Lie algebras are called the classical (Lie) algebras. They are isomorphic to the matrix 
algebras 



Ar^sl(r+1), 5 r ^so(2r + l), C r sp(r), D r ^ so(2r). (2.2) 
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A r series. Let A = (A ss i) be r x r Cartan matrix for the Lie algebra A r = sl(r + 1) , r > 1 
The Cartan matrices for A r -series have the following form 

/ 2 -1 ... \ 
-1 2 -1 ... 
-1 2 ... 



(2.3) 



... 2 -1 
V ... -1 2 ) 

This matrix is described graphically by the Dynkin diagram pictured on Fig. 1. 



1 



r - 1 r 



Fig. 1. Dynkin diagram for A r Lie algebra 

For s s' , A ss i = — 1 if nodes s and s' are connected by a line on the diagram and A ss > = 
otherwise. Using the relation for the inverse matrix A' 1 = (A ss ) (see Sect. 7. 5 in |14| ) 

A ss ' = — |— j- min(s, s') [r + 1 - max(s, s')] (2.4) 

we may rewrite f|l .4[) as follows 

n s = s(r + l-s), (2.5) 

s = 1, . . . , r . 

i? r and C r series. For B r -series we have the following Cartan matrices 



(2.6) 
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while for C r -series the Cartan matrices read as follows 
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Dynkin diagrams for these cases 


are pictured on Fig 


2. 
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Fig. 2. Dynkin diagrams for B r and C r Lie algebras, respectively 
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In these cases we have the following formulas for inverse matricies A 1 = (A ss> ) 



min(s, s') for s 7^ r, 
ns' for s = r, 



A" 



and relation (jl.4p takes the form 



s(2r + 1 — s) for s ^ r, 



for s = r; 



for £? r and C r series, respectively, s = 1, . . . , r . 

D r series. For D r -series the Cartan matrices read 



( 2 -1 
-1 2 -1 
0-12 



min(s, s') for s' ^ r, 
is for s' = r 



n s = s(2r — s), 



\ 





...2 -1 -1 
... - 1 2 
\ 0...-1 2/ 



We have the following Dynkin diagram for this case (Fig. 3): 



(21 



(2.9) 



(2.10) 
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r - 

r - 1 



Fig. 3. Dynkin diagram for D r Lie algebra 
and formula for the inverse matrix A' 1 = (A ss ') [H]: 



A s 



The relation (jl.4p in this case reads 



min(s, s') 


for s, s' ^ {r 


,r-l}, 






for s ^ {r, r 


- 1}, s' G {r,r 


-1}, 




for s £ {r, r 


- 1}, a' i {r,r 


-1}, 


4' 


for s = s' = 


r or s = s' = r 


-1, 




for s = r, s' 


= r — 1 or vice 


versa. 


reads 








J s(2r 


— 1 — s) for 


s £ {r,r - 1}, 






- 1) for 


s £ {r, r — 1}, 





(2.11) 



(2.12) 



s = l,...,r. 

For the simple Lie algebras of type A r , D r , all roots have the same length, and any two nodes 
of Dynkin diagram are connected by at most one line. In the other cases there are roots of two 
different lengths, the length of the long roots being ^/2 times the length of the short roots for B r , 
C T , respectively. 
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3 Computing of fluxbrane polynomials 



At the moment the problem of finding all coefficients Ps in polynomials H s analytically is looking 
as too complicated one. Thereby it is essential to create a program, calculating all the coefficients 
required. 

(k) 

The algorithm of finding coefficients P s is the following one: 

• we substitute the polynomials H s into the set of differential equations (jl.ip and reduce the 

(k) 

differential equations to a set of algebraic equations for P s (expanding the equations into 
degrees of the variable z ) ; 

• the derived system of algebraic equations is solved in terms of the first coefficients P^ = P s ■ 

We choose the symbolic computational system Maple v. 11. 01. for implementation of the dis- 
cussed algorithm. The standard Maple packages LinearAlgebra and PolynomialTools are used 
for working with matrices and polynomials appropriately. So let us start in the following way: 

> with (LinearAlgebra) : 

> with (PolynomialTools) : 

At the next step we need to specify the dimension of the Cartan matrix: 
>S:=3: 

Thus, the dimension of the Cartan matrix is 3. This variable ( S ) also determines the twice 
dual Weyl vector dimension and the number of the differential equations. A variable for the Lie 
algebra is entered similarly. 

>algn : =bn : 

Consequently, by default the program calculates the fluxbrane polynomials for the Lie algebra 
B n with the Cartan matrix, the size of which is 3 . We use the standard Maple function Matrix 
to declare the Cartan matrix. 

>A:=Matrix(S,S) : 

It is more convenient to fill the Cartan matrix with the help of the separate procedure, which 
constructs matrix elements depending on the Lie algebras. Let us consider the procedure. There 
are three callable variables in it: a Lie algebra, the matrix size and the Cartan matrix itself. 

>AlgLie : =proc (algn , S , Cart A : =Matrix (S , S) ) 

The matrix elements are constructed in compliance with Dynkin diagrams. Here we consider 
the construction of the matrix elements for the simple Lie algebras A n , B n , C n and D n . The 
following local variables are essential 

local i,mu,nu; 

i := 0; mu := 0; nu:=0; 
mu := S-l; nu:=S-2; 
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where i is an iteration variable and mu,nu are variables for subscripting in B n , C n and D n 
algebras. By default all elements of the matrix are zeros. Thereby it is necessary to fill only the 
elements different from null according to the Dynkin diagrams. Initially we fill in the main diagonal 
of the matrix, because it is identical to the Cartan matrix of any algebra. 

for i to S do CartA [i, i] := 2 end do; 

Then, using the conditional operator, we consequently take under consideration the condition 
on the matrix size and the fact of belonging to the Lie algebra. We begin with the A n algebra. 

if (S>=1) and algn=an then end if; 

The elements of the secondary diagonals for A n algebra equal —1 according to the Dynkin 
diagram. Firstly, we fill in the upper secondary diagonal, and then it is mirrored one. These actions 
are performed down inside the previous operator. 

for i to S-l do CartA[i, i+1] := -1 end do; 

for i to S-l do CartA [i+1, i] := CartA[i, i+1] end do; 

The similar actions are run with the conditional operator for the B n , in agreement with the 
Dynkin diagram for this algebra. 

if (S>=3) and algn=bn then 

for i from 1 to (S - 1 ) do CartA[i+l,i] :=-l end do; 

for i from 1 to (S - 2) do CartA [i , i+1] : =CartA [i+1 , i] end do; 
end if ; 

But under the Dynkin diagrams for the Lie algebras B n and C n , certain elements are different 
from —1 for the secondary diagonals. So, let us supplement the preceding lines with the following 
matrix element 

CartA [mu,S] :=-2; 

In much the same way we have for the Lie algebra C n 

if (S>=2) and algn=cn then 

for i to S-l do CartA [i, i+1] := -1 end do; 

for i to S-2 do CartA [i+1, i] := CartA [i, i+1] end do; 

CartA [S, mu] := -2; 
end if ; 

For the Lie algebra D n certain non-diagonal elements of the Cartan matrix differ from zero , 
while some elements in the secondary diagonals are equal to zero . According to the Dynkin 
diagram for this algebra the elements of the Cartan matrix are defined in the following way 

if (S>=4) and algn=dn then 

for i to S-2 do CartA [i, i+1] := -1 end do; 

for i to S-2 do CartA [i+1, i] := CartA [i, i+1] end do; 

CartA [S, nu] := CartA [mu, nu] ; 

CartA [nu, S] := CartA [S, nu] ; 
end if ; 
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We return the variable CartA and do not forget to set closing end at the end of the procedure. 
Thus, the Cartan matrix is filled, depending on the chosen algebra, by a call of this procedure with 
the appropriate parameters. 

>AlgLie(algn, S, A) ; 



A :-- 



-1 
2 
-1 





-2 
2 



Further we declare the twice dual Weyl vector by means of the standard Maple procedure 
Vector. 



>n := Vector [row] (1 



S) 



The elements of the inverse Cartan matrix are necessary for calculating the twice dual Weyl 
vector's components. The matrix inversion is done by means of the standard Maple procedure 
Matrixlnverse. 

>A1 := Matrixlnverse (A) ; 



Al := 



1 
2 



We use the standard Maple procedure add for calculating the twice dual Weyl vector's compo- 
nents. 



>for i to S do n[i] := 2*add(Al[i, j] , j = 1 . . S) end do: 

(k) 

The coefficients P s are represented by a matrix. The number of rows of this matrix is 
the number of the differential equations (that is S), and the number of columns is the maximal 
component of the twice dual Weyl vector. But the twice dual Weyl vector was set by means of the 
procedure Vector and the vector must be converted into the list to find the maximal component of 
it using the standard Maple procedure "max" .This action is performed by means of the standard 
Maple procedure convert. 

>maxel := max ( convert (n, list) [] ) : 



Now the matrix of the coefficients P s can be declared. 
>P := array (1 .. S, 1 .. maxel) : 

Let us declare the matrix of the polynomials by means of the procedure Vector. 

>H := Vector [row] (1 . . S) : 

Each element of this matrix is defined according to the hypothesis in the following way: 
>for i to S do H[i] := l+add(P[i, k]*z~k, k = 1 .. n[i]) end do: 
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It's necessary to convert the elements of the matrices H and A into the indexed variables for 
correct calculations. 

>for i to S do for j to S do a[i, j] := A[i, j] end do end do: 
>for i to S do h[i] := H[i] end do: 

Let us enter one more indexed variable o^ v for convenience. 

>for i to S do for v to S do c[i, v] := h[v]~(-a[i, v] ) end do end do: 

We represent the set of equations as a matrix by the use of the procedure Vector. Now the 
system of the differential equations can be defined using the standard Maple procedure dif f . 

>equal := Vector [row] (1 .. S) : 
>for i to S do 

equal[i]:= dif f (z* (dif f (H [i] ,z) ) /H [i] ,z) -P [i , 1] * (product (c [i ,m] , m = 1..S)) 
end do: 

The procedure product is the standard Maple procedure for a product. Further we enter 
two more matrices for simplified equations and numerators of these equations using the procedure 
Vector. 

> simequal := Vector [row] (1 .. S) : 

> newequal := Vector [row] (1 .. S) : 

The first of these matrices is filled, collecting by degrees each of the equations by means of 
the standard procedures. The elements of the second matrix are turned out by selection of the 
numerators from the simplified equations. 

> for i to S do simequal[i] := simplif y(combine(value(equal [i] ) , power)) end do: 

> for i to S do newequal [i] := numer (simequal [i] ) end do: 

(k) 

It's necessary to find out the degrees of the numerators to collect the coefficients P s at 
various degrees of the variable z . So let us describe a matrix, which elements are the degrees of 
the variable z . The degrees are calculated by the standard Maple procedure degree. 

> maxcoeff := Vector [row] (1 .. S) : 

> for i to S do maxcoeff [i] := degree (newequal [i] , z) end do: 

We define a two-dimensional table (the standard Maple structure of data) for the system of 
algebraic equations and fill the table's elements in the following way: 

> coefflist := table(): 

> for i to S do 

for c from to maxcoeff [i] do 

coefflist [i, c] := coeff (newequal [i] , z, c) =0 
end do 
end do : 
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The system should be converted into the list to solve it by means of the standard Maple 
procedure "solve" .This action allows us to apply the function solve. 



> Sys := convert (coeff list , list): 

> sol := solve(Sys) : 

But the form of the answer is inconvenient. Thus substituting the answer into the polynomial 
form, we get 

> trans := {seq(seq(P [i , j] = P[i,j], i = 1..S), j = 1 . .maxel)}: 

> sol := simplif y(map2(subs , trans, sol)): 

> PI := map2(subs, sol, evalm(P)): 

> for i to S do H[i] := l+add(Pl [i ,k] *z~k, k = l..n[i]) end do; 

Hi := 1 + Px,iz + \Pi,xP2,xz 2 + lp 1;1 P 2il P 3il z 3 + -Lp^P^pf ^ + -Lp^pf^pf ^ + 

1 p2 p2 p2 z 6 

129600 1,1 2,1 3,1 ' 

H 2 := 1 + P 2 ,1Z + (^1,1^2,1 + ^P2,lP3,l)z 2 + (^2,1^3,1 + lPl,lP2,lP 3 ,l)z 3 + 
(^ P ll P ll + ^1,1^1,1^3,1 + ^Pl.l^l^l)^ + ^Pl,lP2,lPllZ 5 + 

W^' 1 ^ 1 ' 1 + oTm^ 1 ^ 21 ^ 1 + ^^^^^l^^ 1 ^^ + 32400 ^^i^ 1 ^* 

( - Pi lP 3 P 4 + —L- PhP% X Ph )Z 8 + P 2 P 3 P 4 Z 9 + PhP^Pt Z 10 , 

v 518400 ' 2,1 3l1 259200 1,1 2>1 3 ' 1J 4665600 1,1 2,1 3>1 466560000 1,1 2l1 3,1 
1 



-Pi+pLpLz 6 . 



129600 

It should be noted that throughout the program we use a slightly different notation for the 
first coefficients, i.e. 

P s ,i = P s . (3.1) 

4 Examples of polynomials 

Here we present certain examples of polynomials corresponding to the Lie algebras ^1,^2,^3, 
P3 , C 2 and D4 . 

4.1 ^-polynomials, r = 1,2,3. 

A± -case. The simplest example occurs in the case of the Lie algebra A\ = sl(2) . We get pQ 

H 1 (z) = l + P 1 z. (4.1) 
A 2 -case. For the Lie algebra A 2 = sl(3) with the Cartan matrix 

(A ss ,) = ( _ 2 X ~l ) (4.2) 
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we have [T] n\ = n 2 = 2 and 



H 1 = l + P l z+-P 1 P 2 z 2 , (4.3) 
H 2 = 1 + P 2 z + \PlP2z 2 , (4.4) 



8 = 1,2. 

A 3 -case. The polynomials for the A 3 -case read as follows 



H 1 = l + P 1 z + \P1P2z 2 + ^-P^Psz 3 , (4.5) 
4 36 

F 2 = 1 + P 2 z + (ipxP 2 + \P2Pz) z 2 + ^PiP 2 P 3 z 3 (4.6) 

+—P1P0P3Z 4 , 
144 2 3 ' 

H 3 = 1 + P 3 z + ^P 2 P 3 2 2 + ^^i^Ps^ 3 - (4.7) 
4 3d 



4.2 S3 -polynomials 

For the Lie algebra P 3 we get the following polynomials 

H x = 1 + P x z + -PiP 2 z 2 + — PiP 2 P 3 z 3 + —PiP 2 P%z A + — ^PiP 2 2 PoV (4.8) 
4 18 144 3 3600 2 3 v 7 

+ — - — P?P?P?z 6 , 
129600 1 2 3 ' 

H 2 = l + P 2 z+ (IpyFt + ^3)^ + (^P 2 P| + ^iP 2 P 3 )^ 3 + (^f^f (4-9) 



+— P1P0P3 + — PiP 2 p, 2N )z 4 + — PiP 9 2 Po 2 z 5 + C — !— pip^p 2 + — !— pPpIp 2 

72 1 2 3 16 3 / 600 2 3 V1600 2 3 5184 1 2 3 



1U / UUU \ 1UUU Old 

—P^PDz 6 + f -^—P x PlP* + _^p2p3 2\ 7 / _J_pp3p4 

2592 2 37 V 16200 2 3 32400 1 2 3 J V518400 2 3 

+_^p2p3p|\ 8 + P 2 P 3 P 3 4 z 9 + P?PiPh w , 

259200 1 2 3 / 4665600 1 2 3 466560000 123 

H 3 = l + P 3 z + \p 2 P 3 z 2 + (-^iP 2 P 3 + ^P 3 2 ) z 3 + I Lp 1 P 2 P 2 z 4 (4.10) 

+— *— PiP 2 2 P, 2 ,z 5 + — - — PiP 2 2 P|z 6 . 
3600 2 3 129600 2 3 



4.3 C2 -polynomials 

For the Lie algebra C 2 = so(5) with the Cartan matrix 



-2 2 

we get from (jl.4p rai = 3 and n 2 = 4 . 
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(^')=( I o 1 ! (4-11) 



For C 2 -polynomials we obtain in agreement with [TT] 

Hx = 1 + P x z + \P1P2z 2 + ^PfP2z\ (4.12) 
4 35 

H 2 = l + P 2 z + \P1P2z 2 + \pfP2Z Z + ^-PlPlz\ (4.13) 

4.4 /^-polynomials 

For the Lie algebra D4 we find the following set of polynomials 



H — —P1P0P3P4Z 1 
finn 1 2 3 4 



+ P lZ + ^PP^ 2 + (^1^3 + ^PiP 2 P 4 )^ 3 + -^P^PsPaz 4 (4.14) 

+ -I—P1P0P3P4Z 5 + P?P?P 3 P 4 z 6 , 

3600 2 129600 1 2 

ff 2 = 1 + P 2 z + (-PiP 2 + ^P 2 P 3 + jP2Pa)z 2 + (-P1F2P3 + gPiP 2 Pi (4.15) 

]p 2 p 3 p 4 > 3 + (i^APfPs + ^AP!P4 + -^Pitt + ^PP 2 p 3 p 4 > 4 

? 4 Z 5 + ( P1PIP3P4 + — — PiP 2 2 P. 2 P 4 + -^—P?P?P 3 P A 

V1600 2 5184 2 3 5184 1 2 3 

— !— p 2 p 2 3 p 3 p 4 + — PiP 2 3 p 3 p 4 2 + — : — PiPfp^PiV 

x32400 1 2 32400 2 4 32400 2 3 / 

— - — P?P$P 3 P? + — - — P 2 P 2 3 Po 2 P 4 + — - — P^P^P^z 8 
V518400 1 2 4 518400 1 2 3 518400 2 3 4 / 

+ P 2 PlP 2 Ph 9 + P*P*P*Ph w , 

4665600 1234 46656000 1234 ' 

H 3 = 1 + P 3 z + \p 2 P 3 z 2 + (^PiPiPs + ^PsPi) z * + -^P^PzP^ (4.16) 

+-^—P 1 P 2 2 P 3 P 4 z 5 + — - — PiP 2 2 Po 2 P 4 z 6 , 
3600 1 2 3 4 129600 1234 ' 

H 4 = 1 + P A Z + ^P 2 P 4 Z 2 + (^PlP2Pl + ^^S^ 4 ) ^ + ^A^iV ( 417 ) 

+-^—PiP 2 P 3 P 4 z 5 + — - — PiP?P 3 Ph 6 . 
3600 1 2 3 4 ^ 129600 1 2 3 4 



+— *— PiP?P 3 P?)z 6 + 
5184 2 3 V 



5 Some relations between polynomials 

Let us denote the set of polynomials corresponding to a set of parameters P\ > , P r > as 
following 

H s = H s (z,P 1 ,...,P r ;A), (5.1) 

s = 1, . . . , r , where A = A[Q] is the Cartan matrix corresponding to a (semi)simple Lie algebra 
Q. 
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5.1 C n+ i -polynomials from A2 n +i -ones 

The set of polynomials corresponding to the Lie algebra C n +i may be obtained from the set of 
polynomials corresponding to the Lie algebra ^2n+i according to the following relations 

H s (z, Pi, P n+ i; A[C n+ i]) = H s (z, Pi, P n+ i,P n+ 2 = P n , -F*2n+i = Pi] A[A2 n +i\), (5.2) 

s = l,..., n+1 , i.e. the parameters Pi, P n +i , Pn+2 , P2n+i are identified symmetrically w.r.t. 
P n +i . See Dynkin diagrams on Figs. 1-2. Relation (|5.2p may be verified using the program from 
the Section 3. (For the case n = 1 see formulas from the previous section.) 

5.2 B n -polynomials from /} n+ i-ones 

The set polynomials corresponding to the Lie algebra B n may be obtained from the set of poly- 
nomials corresponding to the Lie algebra D n+ i according to the following relation 

H s (z,Pi,...,P n ;A[B n ]) = H s (z,P u ...,P n ,P n+ i = P n ;A[D n+1 }), (5.3) 

s = 1, . . . ,n, i.e. the parameters P n and P n +i are identified. See Dynkin diagrams on Figs. 2-3. 
Relation (|5.3p may be verified using the program from the Section 3. (For the case n = 3 see 
formulas from the previous section.) 

5.3 Reduction formulas 

Here we denote the Cartan matrix as follows: A = Ay , where T is the related Dynkin graph. Let 
i be a node of T . Let us denote by Tj a Dynkin graph (corresponding to a certain semi-simple Lie 
algebra) that is obtained from T by erasing all lines that have endpoints at i . It may be verified 
(e.g. by using the program) that the following reduction formulae are valid 

H s {z,Pi,...,Pi = 0, ...,P r ;A r ) = H s (z,Pi,...,Pi = 0, ...,P r ;A ri ), (5.4) 
s = 1, . . . , r . Moreover, 

Hi{z, Pi, Pi = 0, .., P r ; A r ) = 1. (5.5) 

This means that setting P{ = we reduce the set of polynomials by replacing the the Cartan 
matrix A-p by the Cartan matrix A^ i ■ In this case the polynomial H \ = 1 corresponds to Ai - 
subalgebra (depicted by the node i ) and the parameter Pi = . 3 . 

As an example of reduction formulas we present the following relations 

H s (z,P 1 ,...,P n ,P n+ i=0;A[g])=H s (z,P 1 ,...,P n ;A[A n ]), (5.6) 

s = l,...,n, for Q = A n+ i, B n+ i, C n +i, D n+ i with appropriate restrictions on n (see (|2.ip ). 
In writing relation (|5.4p we use the numbering of nodes in agreement with the Dynkin diagrams 
depicted on Figs. 1-3. 

The reduction formulas (15. 5p for A§ -polynomials with P3 = are depicted on Fig. 4. The re- 
duced polynomials are coinciding with those corresponding to semisimple Lie algebra ^2 © ^-1 © A2 

• • • • • 

1 2 3 4 5 

Fig. 4. Dynkin diagram for semisimple Lie algebra A2 © Ai © A2 describing the set of 

A5 -polynomials with P3 = 

3 The analytical proof of the relations (|5.2j) - (|5.4[) will be given in a separate publication 
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6 Conclusions 



Here we have presented a description of computational program (written in Maple) for calculation 
of fluxbrane polynomials related to classical simple Lie algebras. (Generalization to semisimple Lie 
algebras is a straightforward one.) This program gives by product a verification of the conjecture 
suggested previously in pp. The polynomials considered above define special solutions to open Toda 
chain equations corresponding to simple Lie algebras that may be of interest for certain applications 
of Toda chains. 

We have also considered (without proof) certain relations between polynomials, e.g. so-called 
reduction formulas. These relations tells us that the most important is the calculation of D„ - 
polynomials, since all other polynomials (e.g. A n -,B n - and C n -ones) may be obtained from 
D n -series of polynomials by using certain reduction formulas. 

A calculation of polynomials corresponding to exceptional Lie algebras (i.e. G2 , i*4 , E$ , Ej 
and Eq ) will be considered in a separate publications. (The Gi -polynomials were obtained earlier 
in [nj.) 
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